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Abstract  R/V  Lance  serendipitously  encountered  an  energetic  wave  event  around  77°N,  26°E  on  2  May  2010. 
Onboard  GPS  records,  interpreted  as  the  surface  wave  signal,  show  the  largest  waves  recorded  in  the  Arctic 
region  with  ice  cover.  Comparing  the  measurements  with  a  spiectral  wave  model  indicated  three  phases  of 
interaction:  (1)  wave  blcxking  by  ice,  (2)  strong  attenuation  of  wave  energy  and  fracturing  of  ice  by  wave  forcing, 
and  (3)  uninhibited  propagation  of  the  fjeak  waves  and  an  extension  of  allowed  waves  to  higher  frequencies 
(above  the  peak).  Wave  properties  during  fracturing  of  ice  cover  indicated  inaeased  groupiness.  Wave-ice 
interaction  presented  binary  behavior  there  was  zero  transmission  in  unbroken  ice  and  total  transmission  in 
fractured  ice.  The  fractured  ice  front  traveled  at  some  fraction  of  the  wave  group  speed.  Findings  do  not  motivate 
new  dissipation  schemes  for  wave  rrxjdels,  though  they  do  indicate  the  need  for  two-way,  wave-ice  coupling. 


1.  Introduction 

The  marginal  ice  zone  (MIZ)  is  the  transition  zone  between  continuous  ice  cover  (so-called  pack  ice)  and  the 
open  oceans.  The  MIZ  is  highly  dynamic  with  coupled  interactions  between  the  ice  floes  and  open  ocean 
processes.  Interactions  between  wind-generated  ocean  surface  waves  and  sea  ice  have  been  a  focus  of 
research  for  some  time  (e.g.,  see  the  reviews  by  Wadhams  [1981],  Sgu/reef  al.  [1995],  and  Squire  [2(X)7]).The 
majority  of  these  studies  have  been  theoretical  in  nature,  based  on  measurements  in  the  laboratory,  or  based 
on  remote  sensing  observations.  To  complement  these  studies,  more  analysis  of  past  in  situ  measurements 
in  the  MIZ  is  beneficial,  particulariy  for  planning  future  field  studies. 

Long  waves  may  propagate  relatively  unhindered  through  ice.  It  is  not  uncommon  for  vessels  to  experience 
swell  events,  even  deep  into  ice  fields,  but  very  few  of  these  events  have  been  reported  over  the  years 
[e.g.,  Robin,  ^963;  Liu  and  Mollo-Christensen,  ‘[988;Asplinet  al.,  201 2],  and  none  of  these  reprart  large  waves 
(e.g,  greater  than  3  m).  The  scarcity  of  in  situ  wave  measurements  in  the  MIZ  is  pierhaps  unsurprising 
considering  the  pro|3ortion  of  global  ship  traffic  that  occurs  in  these  areas.  Very  large  waves  routinely  occur  in 
the  Antarctic  due  to  the  nearly  unlimited  fetch,  yet  measurements  of  large  waves  in  the  MIZ  are  nonetheless 
rare.  Very  recently.  Dob/e  and  B;d/of  [2013]  (henceforth  DB13)  repxjited  buoy  measurements  of  very  large 
waves  in  the  Antarctic  MIZ.  Other  much  less  energetic  (though  still  of  interest)  measurements  have  also 
been  reported  [e.g..  Fox  and  Haskell,  2001;  Meylan  et  al.,  2014].  The  Arctic  fetch  (i.e.,  the  spatial  potential  for 
wave  development)  is  limited  by  ice  cover  (see  the  inset  in  Figure  S5  of  the  supporting  information),  so 
measurement  of  large  waves  is  rarer  still  [Thomson  and  Rogers,  2014].  Several  remarkable  studies  have  been 
made  in  less  energetic  wave  environments  [e.g.,  Hunkins,  1962;  Liu  et  al.,  1991;  Marko,  2003]. 

Mariners  have  long  known  that  wave  energy  diminishes  as  they  navigate  into  ice,  and  many  have  witnessed 
swell  propjagate  deep  (distances  of  2  to  4  orders  of  magnitude  greater  than  the  characteristic  wavelength) 
into  the  ice.  Research  has  focused  on  several  effects  of  ice  cover  on  incoming  waves:  change  in  dispersion 
relation  [e.g.,  Hunkins,  1 962;  Liu  ef  o/,  1 991 ;  Fox  and  Haskell,  2001  ]  and  other  radiative  effects  such  as  reflection, 
transmission  [GoMshteh  and  Marchenko,  1989;  Fox  and  Squire,  1994],  scattering,  refraction,  and  attenuation 
[Squire  et  al.,  1 995;  Perrie  and  Hu,  1 996;  Marchenko  artd  Voliak,  1 997;  Liu  and  Mollo-Oiristensen,  1 988;  Squire  and 
WilSams,  2008].  The  wave  effects  on  ice  are  mostly  mechanical  in  nature:  flexing  and  fracturing  of  continuous  ice 
and  floes,  calving  of  ice  edges,  convergence  or  divergence  of  ice  fields,  and  forcing  collisions  between  floes 
[Wadhams,  1981;  Squire,  2007,  and  references  within].  Ice  is  often  heterogeneous  in  nature  and  varied  in  form 
which  greatly  impacts  wave-ice  interaction  [e.g.,  Campbell  et  al.,  2014]. 
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Of  sp>ecial  interest  is  the  fracturing  and  convergence  of  ice  floes  under  the  influence  of  swell  because  this 
represents  a  possible  feedback  loop  between  air,  sea,  and  ice  [Asplin  et  al.,20M;Thomson  and  Rogers,  2014; 
Asplin  et  al,  201 4].  On  the  macroscale  (e.g.,  seasonal),  action  of  wind  on  the  open  sea  generates  waves; 
these  waves  contribute  to  the  fracturing  of  ice  through  mechanical  strain  (i.e.,  creeping  and  bending).  By 
increasing  the  surface  area  of  individual  ice  floes,  this  fracturing  accelerates  melting.  The  newly  ice-free 
regions  increase  the  area  of  open  sea  (i.e.,  fetch),  in  turn  allowing  larger  waves  to  form  and  break  up  more 
ice.  it  is  thought  that  this  process  will  become  more  important  as  the  warming  climate  reduces  the  extent 
of  perennial  Arctic  ice  cover  [Squire.  2007;  Squire  et  al,.  2009;  Khon  et  al.,  2014).  On  the  miaoscale  (time  scale 
of  hours),  we  propose  that  large-wave  events  fracture  ice,  creating  an  ice  cover  that  permits  the  passage  of  a 
larger  fraction  of  the  wave  energy  (as  documented  by  OBI  3),  leading  to  more  wave  energy  further  into  the 
ice;  arrd  further  fracturing  of  the  ice,  until  an  equilibrium  (between  wave  strain  and  ice  strength)  is  reached  or 
the  event  ends,  and  the  event-induced  migration  of  the  MIZ  toward  the  interior  is  halted. 

Traditionally,  sea  ice  has  been  treated  simplistically  in  the  spectral  models  that  are  used  for  routine  hirxicasting 
and  forecasting  of  surface  waves,  such  that  the  MIZ  is  treated  as  either  op>en  water  or  imprermeable  (as  land) 
based  on  an  ice  concentration  criterion  [Tuomi  et  al,  201 1].  This  was  improved  somewhat  by  Tolman  [2003] 
using  a  cell-wise  partial  blocking  of  energy  flux,  scaled  by  ice  concentration.  Wave-ice  interaction  remains  an 
area  of  wave  modeling  which  is  ripe  for  progress.  More  recently,  efforts  have  been  made  to  incorporate  the 
effect  of  ice  into  the  physics  of  these  wave  models  [DB1 3;  Rogers  and  Orzech,  201 3;  Rogers  and  Zieger,  2014; 
Tolman  and  the  WAVEWATCH  III*  Development  Group,  2014].  Wave  energy  can  be  either  scattered  and  reflected 
(conservative  processes)  or  dissipated  (a  nonconservative  process).  Either  causes  a  reduction  of  wave  energy 
along  the  main  axis  of  propagation.  With  respect  to  dissipation,  there  are  a  few  theoretical  paradigms,  any  of 
which  may  be  valid  under  certain  drcumstances  owing  to  the  varied  nature  of  ice  floes.  Two  that  have  been 
proposed  are  dissipation  through  turbulence  generated  beneath  ice  [Liuetal.,  1991]  and  dissipation  through  an 
effective  viscosity  of  the  ice  cover  [Keller,  1 998;  Newyear  and  Martin,  1 999;  Wang  and  Shen,  2010].  In  the  case  of 
the  implementation  in  the  WAVEWATCH  III*  model  [Rogers  and  Orzech,  2013;  Rogers  and  Zieger,  2014;  Tolman 
and  the  WAVEWATCH  III*  Development  Group,  201 4],  the  new  formulations  remained  uncalibrated  except  in  the 
most  gross  sense;  and  therefore,  field  studies  with  detailed  wave  measurements  (such  as  the  one  desaibed 
herein)  are  seen  as  an  opportunity  to  perform  such  calibrations.  The  data  set  of  the  present  study,  though 
not  ideally  suited  for  such  calibration  work,  does  benefit  our  understanding  of  wave-ice  interaction  in  terms  of 
model  requirements  (see  section  4). 

Some  of  the  aforementioned  reports  arose  from  data  of  opportunity,  that  is  to  say  the  field  experiments 
were  not  designed  for  wave  measurements,  and  the  same  may  be  said  for  our  study.  To  help  fill  the  gap  of 
in  situ  observations  of  wave-ice  interaction,  we  present  the  largest  waves  measured  in  the  Arctic  under 
significant  ice  cover  (section  3)  and,  during  the  breakup  phase  (section  32),  closely  examine  wave  properties 
and  infer  ice  quality.  These  measurements  are  compared  with  results  from  a  spectral  wave  model  (section  3), 
and  the  requirements  for  accurate  wave  prediction  are  discussed  (section  4).  In  addition,  the  possible  role  of 
nonlinear  wave  focusing  during  ice  breakup  is  discussed  (section  4.1). 

2.  Methods 

The  R/V  Lance  is  a  former  fishing  vessel,  converted  for  research  and  expedition  and  operated  by  the 
Norwegian  Polar  Institute.  This  particular  cruise  was  organized  for  the  field  works  section  of  the  University 
Center  in  Svalbard's  course  concerning  drifting  ice  in  the  Barents  Sea.  The  R/V  Lance  is  60.8  m  in  length  by 
12.6  m  in  breadth.  The  location,  including  vertical  position,  of  the  vessel  was  recorded  by  an  onboard 
Seapath  200  GPS  system  at  1  Hz.Theaccuracyof  the  vertical  position  was  greater  ofO.OSm  or  5%  (defined 
in  the  Seapath  200  technical  documentation).  Using  this  record,  the  ship  velocity  was  determined  from 
tracking  the  location  and  the  sea  surface  elevation,  q,  from  detrending  the  vertical  position.  Treating  the 
vessel  as  a  surface  following  buoy  whose  dimensions  act  as  a  low-pass  filter  to  the  signal  from  the 
ambient  sea  surface  spectrum,  1-D  frequency  spectra  were  calculated  from  1  h  records  of;/.  In  theory,  if 
the  vessel  is  stationary,  waves  with  lengths  at  least  twice  its  dimensions  may  be  accurately  resolved 
assuming  a  perfect  response.  This  implies  a  minimum  wavelength  (maximum  frequency  assuming  linear 
disprersion  relation)  of  121 .6  m  (0.1 13  Hz)  for  waves  approaching  head  on  and  25.2  m  (0.249  Hz)  for  waves 
approaching  the  broadside. 
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Since  R/V  Lance  is  most  effective  in  measuring  relatively  low  frequency  wave  energy,  we  will  compare  an 
equivalent  low-frequency  height  by  integrating  the  1-D  wave  spectrum,  S(f),  in  analogy  to  significant  wave 
height  but  with  an  upper  frequency  bound 


Later  in  the  p»af>er,  we  discuss  wave  groupiness  (section  4.1 ).  We  calculated  two  associated  metrics:  the  so-called 
Benjamin-Feir  Index,  BFI  [Janssen,  2(X)31  and  the  relative  sp>ectral  width  parameter,  v  [Babanin  and  Sobviev,  1 998]: 


BFI  =  2y/23rkpn^m~^'^ 


where  fp  is  the  frequency  at  the  peak  (i.e.,  most  energetic  region)  of  the  spectrum  (aka  peak  frequency)  and  kp 
is  determined  by  the  deep  water  linear  dispersion  relation,  fp  =  ^gkpjhc.  The  usual  zeroth-order  spectral 
moment  is  mo,  and  n  is  a  modified  spectral  moment  as  follows: 


We  do  not  compare  these  parameters  with  a  model,  hence,  the  spectral  bounds  have  been  relaxed  to  include 
higher-frequency  waves.  Wave  measurements  from  a  vessel  underway  present  a  problem  for  interpretation 
of  the  signal.  This  is  because  the  forward  speed,  with  no  changes  in  velocity  or  heading,  enters  into  the 
dispersion  relation  as  an  extra  term: 


This  is  the  weli-known  Doppler  shift  (e.g.,  Hanson  et  a!.,  1 997]  which  changes  the  frequency,  measured  in  the 
ship  frame  of  reference,  according  to  the  (1 )  relative  angle  between  the  wave  propagation  direction  and  the 
ship  heading,  a,  and  the  (2)  product  of  the  magnitudes  of  the  ship  velocity,  V,  and  wave  number,  k  (derived  from 
the  intrinsic  frequency).  To  avoid  this  issue,  analysis  was  confined  to  a  period  when  the  ship  speed  was  below 
3  m/s  (note  that  the  frequencies  considered  in  have  phase  speeds  >  1 3  m/s). 

A  nested  run  of  SWAN  ('Simulating  WAves  Nearshore")  [Booij  et  a!.,  1999]  was  performed  as  part  of  the 
analysis  of  this  wave  event.  The  outer  computational  grid  was  defined  from  [  24°E,  66°N]  to  [SS'E,  79°N] 
with  ~1/3°  resolution  and  the  inner  computational  grid  from  [15°E,  74°N]  to  [30°E,  78°N]  with  ~1/10°  resolution. 
The  model  topography  was  from  ET0P02  [National  Geophysical  Data  Center,  2006L  arvd  the  10  m  wind  vectors 
were  obtained  from  archives  of  the  U5.  Nav/s  operational  analyses  (Navy  Operational  Global  Atmospheric 
Prediction  System)  [Hogan  and  Rosmond,  1991].  SWAN  was  run  without  any  representation  of  the  icci  though 
operational  analyses  from  the  Navy's  Polar  Ice  Prediction  System  (PIPS)  [Van  Woertetal.,  2004]  were  used 
to  give  a  large-scale,  spatial  overview  of  the  concentration  of  sea  ice.  PIPS  was  updated  once  per  day  at 
noon  UTC.  For  local,  qualitative  inferences  of  ice  conditions,  we  utilized  the  first-hand  account  of  author 
A.M.  and  bow-oriented  photographs  captured  hourly. 

3.  Results 

3.1.  Phase  1 — ^Wave  Blocking  by  Solid  Ice 

Phase  1  is  marked  in  red  in  Figure  1  and  spans  the  time  period  from  2  May  2010  15:00-20:15  LTFC.  On  1  May 
2010,  a  low-pressure  system  over  northern  Europe  moved  northeast  and  strengthened  over  the  Barents  Sea 
on  2  May  2010.  During  this  time,  field  work  required  R/V  Lance  to  moor  to  pack  ice  about  75  km  north 
northeast  of  Hopen  Island  (Svalbard,  Nonway).  Visual  ice  observations  indicated  continuous,  p»ack  ice  with 
relatively  flat  surfaces  covered  by  snow.  Ridge  sail  heights  did  not  exceed  1  m,  and  the  ice  thickness  was 
typical  for  this  area  at  this  time  of  year  (~05-0.6  m).  The  rigid  structure  of  the  ice  can  be  seen  in  Figure  1  d. 
Since  there  was  no  ice  representation  in  SWAN,  SWAN  predicted  a  gradual  increase  of  wave  height  from 
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Figure  1 .  (a)  Map  of  southern  Svalbard  in  the  top  left  and  Hopen  Island  in  the  bottom  right.  The  ship  track  is  in  color  corresponding  to  different  phases  of  wave 
interaction  (red:  wave  blocking  by  ice,  yellow:  arrival  of  waves  and  breakup  of  the  ice,  green:  swell  event  at  full  energy  allowing  higher  frequencies  with  time, 
cyan:  the  swell  shadow  of  Hopen  Island,  blue  and  purple:  after  the  swell  shadow).  Black  dots  correspond  In  space  to  the  black  dots  in  the  time  series.  Stars 
(and  colors)  correspond  with  the  spectra  shown  in  Figure  3.  The  background  is  superimposed  with  an  advanced  synthetic  aperture  radar  image  captured  on  1 
May  2010  09:45  38  UTC  on  the  Envisat  satellite  operated  by  European  Space  Agency  and  provided  by  Benjamin  Holt  at  the  Jet  Propulsion  Laboratory,  (b)  Time 
series  of  H|ow  as  measured  by  IWLonce  (black  solid  line)  and  predicted  by  SWAN  (black  dashed  line)  and  the  ship  velocity  (grey  solid  line).  Stars  correspond  with 
the  spectra  shown  in  Figure  3.  (c)  BFI  (black  solid  line)  and  v  (black  dashed  line),  (d  i)  Select  photographs  from  the  ship  show  local  sea  ice  state.  A  color  bar, 
corresponding  to  the  colors  of  ship  track,  is  located  at  the  bottom  of  each  time  series  for  convenience.  The  phases  of  wave  interaction  are  labeled. 


around  1  m  at  -0900  to  4  m  at  ~21:(X),  during  which  time  the  pack  ice  was,  in  reality,  blocking  almost  all 
incoming  wave  energy.  The  first  waves  were  measured  at  ~1 9:40  as  an  isolated  wave  group  (see  Figure  2),  the 
central  wave  height  of  which  was  0.9  m.  Soon  after,  newly  formed  cracks  in  the  ice  were  discovered,  and  it 
was  deemed  too  dangerous  to  continue  field  operations.  It  is  probable  that  wave  energy  was  transmitted  into 
the  ice  as  flexural-gravity  waves,  of  amplitude  too  small  for  the  ship  to  measure  reliably,  but  which  were 
responsible  for  the  formations  of  the  cracks  in  the  ice. 

3.2.  Phase  2 — Ice  Breakup  and  Transition 

Phase  2  is  marked  in  yellow  in  Figures  1  and  2.  At  -20:15,  less  than  an  hour  following  the  isolated  wave  group, 
the  ship  again  encountered  small  waves.  Over  the  course  of  the  next  hour,  the  aacked  ice  gave  way  and  broke 
into  smaller  floes,  and  the  observed  waves,  which  were  initially  under  1  m,  very  quickly  became  4  m  swell. 
Figure  2e  shows  the  time  series  of  tj  and  photographs  of  the  ice  during  the  transition.  The  ice  is  largely 
unchanged  between  Figures  2a  and  2b.  In  Figure  2c,  the  ice  appears  slightly  more  homogenous  while  still 
retaining  larger  floes,  while  in  Figure  2d,  the  ice  is  clearly  broken  into  smaller,  more  uniform,  floes.  The  waves 
during  the  transition  perbd  have  a  grouped  (modulated)  structure  in  Figure  2,  which  manifests  as  elevated  BFI 
and  lower  v  in  Figure  1c.  At  the  end  of  this  1  h  transition  period,  the  Huiw  matched  the  prediction  from  SWAN. 

3.3.  Phase  3 — Continued  Fracturing 

Phase  3  is  marked  in  green  in  Figures  1  and  2.  The  ship  continued  to  measure  large  swell  over  the  next  7  h  as  it 
slowly  steamed  through  the  MIZ  headed  for  protected  waters  behind  Hopen  Island.  Simple  analysis  of 
Figure  1  g  (see  supporting  information)  showed  that  the  peak  wavelength  (-200  m)  was  much  longer  than 
the  typical  floe  size  (-5-10  m).  H|ow  from  SWAN  agreed  nearly  pierfectly  with  the  ship  measurements  during 
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2  May  2010 


Figure  2.  Time  series  of  sea  surface  elevation,  >),  derived  from  GPS  on  RAf  Lance  and  corresponding  photographs  of  the  ocean  surface  before,  during,  and  after 
the  ice  breakup. 


this  period,  implying  that  the  ice  had  little  or  no  effect  on  the  waves  measured  by  the  ship.  The  measured 
sprectra  in  Figure  3,  which  correspond  to  the  stars  in  sp>ace  and  time  in  Figures  2a  and  2b,  respectively,  show 
that  there  was  a  continuous  increase  in  energy  at  frequencies  above  the  SF>ectral  peak.  At  ~04:00  on  3  May 
2010,  the  ship  reached  the  swell  shadow  of  Flopen  Island  and  allowed  for  breakfast  before  steaming  onward, 
at  an  increased  pace,  to  the  south.  The  vessel  again  encountered  large  waves  at  ~08:00  (seen  as  blue  in 
Figure  1 )  as  it  exited  the  swell  shadow  headed  east. 

4.  Discussion 

4.1.  Groupiness 

The  waves  which  were  responsible  for  the  initial  breakup  of  the  pack  ice  were  "groupy,”  meaning  the  wave 
energy  arrived  in  discrete  packets  or  groups  (also  referred  to  as  "sets,"  "beats,"  or  "envelopes").  This  is  evident 
in  the  time  series  of  tj  in  Figure  2  and  is  corroborated  by  spectral  parameters  BFI  and  v  in  Figure  1c.  BFI  is 

the  ratio  of  nonlinearity  to  dispersion, 
and  dispersion  is  related  to  the  spectral 
width.  From  Figure  3,  it  is  evident  that 
the  spectra  were  filtered  (more  on  this 
follows  below),  effectively  decreasing 
the  spectral  width  which  manifest  as 
a  larger  (smaller)  BFI  (v)  near  the 
transition  (phase  2).  This  indicates 
that  the  wave  system  may  have  been 
subject  to  modulation  instability  (Ml) 
[Benjamin  and  Feir,  1967;  see  also  Yuen 
and  Lake,  1 982].  It  was  shown  by  Liu 
and  Mollo-Christensen  [1 988]  that 
monochromatic  waves  more  quickly 
destabilize  into  groups  when  the  sea 
surface  is  under  high  compression 
stress  (i.e,  ice  cover)  than  over  the 


10'. 


Spectral  density 


Figure  3.  The  evolution  of  S(f)  over  the  course  of  5  h  as  the  ship 
encountered  smaller  ice  floes.  The  first  spectra  is  shown  in  dark  blue,  and 
later  spectra  transitioning  to  aqua.  SWAN  spectra  at  this  time,  with  no  ice 
representation,  is  shovm  in  black  for  reference. 
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open  ocean.  If  this  is  true,  the  nonlinear  focusing  of  wave  energy  via  Ml  may  possibly  explain  how  large 
waves  are  able  to  p)enetrate  more  deeply  into  ice  cover.  A  feedback  mechanism  is  possible:  the  compact 
nature  of  the  (ice  covered)  sea  surface  provides  a  medium  which  enhances  the  self  focusing  of  the  waves 
(i.e..  Ml),  these  waves,  in  turn,  increase  the  strain  which  leads  to  the  initial  breaking  of  the  ice,  broken  ice 
allows  the  waves  to  travel  more  deeply  into  covered  ice.  in  this  way,  nonlinear  wave  trains  lead  the  advance 
into  the  unbroken  ice  field.  After  the  initial  passage  of  wave  groups,  ice  no  longer  provides  the  type  of 
stress  which  leads  to  enhanced  self  focusing,  nor  does  it  inhibit  freely  propagating  linear  waves  in  its 
fractured  state.  This  transition  period  is  about  1  h  in  length  for  our  wave  event.  This  represents  a  new  and 
potentially  important  mechanism  for  rapid  ice  breakup. 

This  theory,  however,  must  be  temprered  by  uncertainties  in  the  explanation.  Most  significant  is  the  exact  nature 
of  the  disprersion  relation  (imprortant  for  the  calculation  of  BFI  via  steepness)  in  the  MIZ  In  ice,  the  disprersion 
relation  has  additional  deprendencies  on  ice  proprerties  including  thickness,  eicistic  modulus,  compression  stress, 
and  viscosity.  The  amount  of  compression  stress  assumed  in  Liu  and  Mollo-Christensen  [1988]  may  not  be 
representative  of  typical  sea  ice  conditions  [e.g.,  Hibler,  1 979,  also  see  Timco  and  Weeks,  201 0].  The  detailed 
ice  properties  during  phase  2  in  this  study  are  unknown;  furthermore,  the  exact  nature  of  the  dispersion 
relation  in  ice  is  a  matter  of  contention  for  which  there  is  little  observational  evidence.  Certainly,  much 
more  research  is  needed  to  establish  the  imprortance  of  such  a  mechanism. 

4.2.  Ice  as  a  Low-Pass  Filter 

Figure  3  shows  the  evolution  of  wave  spectra  with  time.  Measured  spjectra  are  shown  in  color  indicating  the  time, 
and  a  spectrum  from  SWAN  is  shown  in  black  as  a  "no  ice"  reference.  The  sp>ectral  density  changes  little  below 
the  spectral  p)eak,  but  above  the  peak,  more  and  more  energy  was  being  detected  as  time  progressed.  This 
observation  is  similar  to  that  of  a  low-|3ass  filter  in  which  the  cutoff  frequency  is  increasing  with  time.  We  suggest 
that  the  ice  field  acted  as  the  filter  [Marko,  2003].  As  evidenced  in  the  photographs  in  figures  2  arxJ  3,  the  ice  field 
was  characterized  by  smaller  floes  as  time  passed.  The  smaller  dimensions  of  the  floes  allowed  higher-frequency 
waves  to  pxopjagate.  Similar  features  may  be  seen  in  the  observations  of  Robin  [1963,  Figures  7  and  8]. 

The  smaller  floe  diameters  may  have  been  a  result  of  wave  action  (temporal  change  in  ice  structure)  or  the 
ship  moving  toward  the  open  water  through  the  MIZ  (spjatial  change  in  ice  structure).  The  satellite  image  in 
Figure  1  a  and  the  ice  concentrations  from  PIPS  (see  Figure  S5  of  the  supporting  information)  suggest  the 
latter  prossibility.  In  all  likelihood,  it  was  a  combination  of  nonstationary  and  nonuniform  ice  contributing  to 
the  varied  ice  characteristics  observed  by  the  moving  ship. 

4.3.  Fractured  Ice  Front 

We  hyprothesize  that  there  was  a  front  seprarating  nearly  solid  prack  ice  and  fractured  ice. This  front,  driven  by 
wave  action,  was  traveling  deep>er  into  the  ice  field.  Eventually,  this  front  crossed  the  p>osition  of  the  ship,  the 
time  at  which  it  first  encountered  the  wave  event.  SWAN  was  used  to  estimate  the  delay  of  swell  event  to  be 
~1 2  h.  PIPS  was  used  to  estimate  the  distance  from  the  position  of  R/V  Lance  to  the  closest  ice  edge  in  the 
direction  of  wave  propagation  to  be  ~200  km.  With  these  two  quantities,  the  sp>eed  of  the  front  is  estimated 
to  have  been  roughly  45  m/s,  which  happens  to  be  about  half  the  wave  group  velocity.  The  fact  that  the 
estimated  speed  of  the  fractured  ice  front  was  some  significant  fraction  of  the  group  sp>eed  means  the  nature 
of  the  icefield  was  changing  (in  time  and  space)  on  scales  similar  to  those  which  characterize  the  peak  waves. 

4.4.  Implications  for  Ice  Implementation  in  Spectral  Wave  Models 

The  results  presented  here  have  implications  for  spectral  wave  models,  although,  being  a  case  study,  the 
generality  of  the  implications  is  not  known.  There  is  limited  evidence  that  large  waves  behave  differently  in 
ice  [Kohout  et  al.,  20 1 4];  as  such,  it  is  unclear  if  this  study  represents  typical  wave-ice  interactions.  Kohout  et  ai. 
[2014]  claim  that  wave  heights  greater  than  3  m  (for  a  certain  wave  period)  are  linearly  attenuated  by  ice, 
whereas  waves  of  heights  less  than  3  m  are  exponentially  (i.e.,  more  strongly)  attenuated.  Our  results  also 
support  less  attenuation  for  larger  waves.  In  fact,  the  close  agreement  with  the  model  (keeping  in  mind  there 
is  no  ice  representation  in  the  model)  implies  there  is  negligible  wave  attenuation  at  the  preak  frequencies 
by  fractured  floes.  This  is  in  agreement  with  results  of  DB1 3.  Our  observations  indicate  wave-ice  interaction 
with  binary  behavior,  where  peak  waves  are  completely  blocked  (i.e.,  100%  dissipation),  then  completely 
p)assed  (i.e,  zero  dissip>ation)  with  a  very  short  transition  time  (0(1  h)).  This  typo  of  behavior  has  important 
implications  for  wave  modeling.  If  the  wave  model  was  "aware”  of  the  location  of  the  solid  ice,  it  could  well 
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predict  the  wave  conditions  by  not  allowing  energy  to  propagate  into  that  region.  We  have  shown  that  it 
may  be  possible  for  the  location  of  the  solid  ice  edge  to  move  at  some  significant  fraction  of  the  wave 
group  speed.  Hence,  the  wave  model  would  not  perform  well  with  information  from  routine  ice  products 
which  are  low  resolution  in  space  (0(1 0  km))  and  sparse  in  time  (updated  once  per  day).  It  is  also  not  clear 
that  ice  concentration  would  necessarily  be  different  for  unbroken  and  broken  ice,  so  perhaps,  a  new  type 
of  ice  product  is  needed.  In  essence,  accurate  wave  prediction  would  require  two-way  coupling  with  a 
high-resolution  model  which  could  resolve  such  processes  as  described  here.  The  skill  of  models  to  forecast 
waves  in  the  Arctic  is  then  crucially  dependent  on  the  accuracy  of  the  ice  forecast,  requiring  resolution  of 
ice  fields  on  time  and  space  scales  of  hours  and  kilometers.  In  reality,  operational  implementation  at  such 
scales  is  not  presently  possible. 

4.5.  Climate  Change? 

Singular  events  cannot  be  attributed  to  long  time  scale  processes  such  as  climate  change.  So  the  question  of 
whether  or  not  this  large  wave  event  in  Arctic  ice  is  a  consequence  of  a  changing  dimate  is  ill  posed.  A  better 
question  is,  will  energetic  wave  events  occur  more  frequently  in  the  Arctic  in  the  future  as  the  climate  warnis? 
Our  answer  is  yes,  a  warmer  dimate  will,  as  one  reviewer  noted,  lead  to  dedining  summer  sea  ice  cover  (arnd 
delayed  freeze-up  in  the  fall).  While  fetch  in  the  Nonwegian  and  Barents  Sea  has  historically  been  available,  the 
extent  of  ice-free  regions  has  recently  reached  historic  levels,  and  completely  new  fetches  have  been  exposed  in 
the  Beaufort,  Chukchi,  and  East  Siberian  Seas.  Furthermore,  there  is  evidence  of  increased  storm  activity  in  the 
Arctic  [Sepp  and  Jaagus,  201 1  ].  All  of  this  effectively  exacerbates  the  fetch-wave,  wave-ice,  ice-fetch  feedback 
loop,  and  highly  energetic  wave  events  will  occur  more  frequently  [e.g,  Dobrynin  etal,  2012;  Khon  eta!.,  2014]. 

A  few  of  the  potential  impacts  of  increased  wave  heights  may  be  rerouting  of  shipping,  damage  to  industrial 
platforms,  and  the  loss  of  land  in  coastal  communities  to  erosion.  As  interest,  presence,  and  research  in  the 
Arctic  intensifies,  these  observations  are  soon  to  be  complemented  by  measurements  of  more  and  more 
energetic  wave  events. 
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5.  Conclusions 

An  energetic  wave  event  under  significant  ice  cover  in  the  Arctic  was  measured  by  R/V  Lance  on  2  and  3  May 
2010  in  the  vicinity  of  Svalbard,  Nonway.  The  results  indicate  three  distinct  phases  of  ice-wave  interaction. 
(Xiring  phase  1,  the  waves  are  completely  blocked  and/or  scattered  by  ice  and  a  spectral  wave  rrrodel  (without 
any  ice  representation)  overpredicts  the  wave  height,  as  expected.  During  phase  2,  the  ice  was  fractured  by 
wave  forcing  and  the  fractured  ice  began  to  allow  wave  energy  to  propagate.  This  transition  period  lasted  about 
an  hour,  during  which  time  the  waves  appeared  mcxJulated  or  groupy.  Mcxlulaticjn  instability,  hypothesized 
to  be  triggered  by  a  change  in  wave  dynamics  in  ice,  is  prcjposed  as  a  possible  mechanism  for  the  initial  breakup 
of  ice.  Firm  conclusions  were  hindered  by  the  lack  of  detailed  ice  properties  and  the  speculative  nature  of 
the  dispersion  relationship  in  ice  covered  seas,  but  the  authors  regard  this  as  a  promising  avenue  for  future 
research.  During  phase  3,  the  wave  model  and  measurements  were  in  excellent  agreement,  implying  peak 
waves  propagated  through  the  ice  uninhibited.  The  agreement  continued  over  the  course  of  7  h  during  which 
time  the  measurements  indicate  an  extension  of  allowed  frequendes  to  those  above  the  peak.  Our  results 
suggest  that  accurate  wave  prediction  would  have  required  coupling  with  an  ice  model  which  resolved  scales  of 
hours  and  kilometers.  Implementation  at  such  scales  will  no  doubt  be  a  future  challenge.  These  are  the  largest 
known  waves  recorded  in  the  Arctic  with  substantial  ice  cover  present,  and  we  exprect  the  measurement  of 
large-wave  events  to  occur  more  frequently  in  the  future  due  to  the  fetch  wave-ice  fetch  feedback  loop. 
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